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1^ ■ Abstract 
m . 

Q"^ . We investigate target mass effects in tfie unpolarized virtual photon structure functions 



F2{x,Q'^,P'^) and F2{x,Q'^ , P'^) in perturbative QCD for the kinematical region ^ ^ Q^, 
where —Q'^{—P'^) is the mass squared of the probe (target) photon and A is the QCD scale param- 
eter. We obtain the Nachtmann moments for the structure functions and then, by inverting the 



^ I moments, we get the expressions in closed form for -Fg (^) Q ^P ) up to the next-to-next-to- leading 



d ' order and for F7(x, Q^, P^) up to the next-to-leading order, both of which include the target mass 



corrections. Numerical analysis exhibits that target mass effects appear at large x and become 
sizable near a;inax(= 1/(1 + ^))' maximal value of the ratio P^/Q^ increases. 
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I. INTRODUCTION 



It is well known that, in e~^e~ collision experiments, the cross section for the two-photon 
processes e"*"e~ —>■ e^e~ +hadrons illustrated in Fig. 1 dominates at high energies over other 
processes such as one-photon annihilation process e~^e~ — 7* — hadrons. Here we consider 
the two-photon processes in the double-tag events, where both the outgoing e"*" and e~ are 
detected. Especially, we investigate the case in which one of the virtual photon is far off-shell 
(large = —q"^), while the other is close to the mass-shell (small = —p"^). This process 
can be viewed as a deep-inelastic scattering off a photon target \^ with mass squared —P^, 
through which we can study the photon structure functions. 

In the case of a real photon target (P^ = 0), unpolarized (spin-averaged) photon structure 
functions F2{x,Q^) and Fl{x,Q'^) were studied first in the parton model p.], and then 
investigated in perturbative QCD (pQCD). In the framework based on the operator product 
expansion (OPE) [s] supplemented by the renormalization (RG) group method, Witten 4] 
obtained the leading order (LO) QCD contributions to F2 and Fl and, shortly after, the 
next-to-leading order (NLO) QCD corrections to FJ were calculated by Bardeen and Buras 

n HQ 

[5|]. The same results were rederived by the QCD improved parton model approach p, 17|. 
The QCD analysis of the polarized photon structure function gl{x,Q'^) for the real photon 
target was performed in the LO sl] and in the NLO 

The structure functions Fiji^x, Q^, P^) and F2{x, Q^, P^) for the case of a virtual photon 
target (P^ ^ 0) were studied in the LO and in the NLO 13] by pQCD. In fact, these 
structure functions were analyzed in the kinematical region, 

A' < P' < , (1.1) 

where A is the QCD scale parameter. The advantage of studying a virtual photon target in 
the kinematical region (11.11) is that we can calculate the whole structure function, its shape 
and magnitude, by the perturbative method. This is contrasted with the case of the real 
photon target where in the NLO there exist nonperturbative pieces. The virtual photon 
structure functions Fn and F2 were also studied by using the DGLAP-type QCD evolution 



equations 
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16| . In the same kinematical region (11. II) . the polarized virtual 



photon 



m the second paper of [4 
photon were analyzed in 



17| and 



structure function g1{x , Q"^ , P"^) was investigated up to the NLO in QCD in Ref. 

I. Moreover, the polarized parton distributions inside the virtual 
18|. Recently the first moment of gJ{x,Q^,P'^) was calculated 



e-(e^) 




FIG. 1: Deep inelastic scattering on a virtual photon in the e"*" e collider experiments. 



up to the next-to-next-to-leading order (NNLO) [19]. For more information on the recent 
theoretical and experimental investigation of unpolarized and polarized photon structure, 
see the review articles {2^. 

In our previous paper j2l|, we have studied the unpolarized virtual photon structure 
functions, F^{x,Q\P^) up to the NNLO and F2{x,Q\P^) up to the NLO, in pQCD 
for the kinematical region fll.ll) . This investigation became possible thanks to the recent 
:hree-loop calculations of the parton-parton as well as photon-parton splitting functions 



22I . I23I . |2J]. There we have considered the logarithmic corrections arising from the QCD 
higher-order effects up to the NNLO, and ignored all the power corrections of the form 
{P'^ IQ'^Y = 2, ■ ■ ■ ) coming either from target mass effects or from higher-twist effects. 

In fact, if the target is a real photon (P^ = 0), there is no need to consider target mass 
corrections. But when the target becomes off-shell, for example, P^ > M^, where M is 
the nucleon mass, and for relatively low values of Q^, contributions suppressed by powers 
of P^/Q^ may become important. Then we need to take into account these target mass 
contributions just like the case of the nucleon structure functions. The consideration of 
target mass effects (TME) is important by another reason. For the virtual photon target, 
the maximal value of the Bjorken variable x is not 1 but 

_ 1 

2^max ~ ) (-'-•2) 

due to the constraint {p + > 0, which is in contrast to the nucleon case where Xmax = 1- 
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The structure functions should vanish at x = Xmax- However, both the QCD NNLO result for 
F^i^x, Q^, P^) and the NLO result for Fl{x, Q^, P^) [21[ show that the predicted graphs do 
not vanish but remains finite at x = Xmax- This flaw is coming from the fact that TME have 



not been taken into account in the analysis 
in the past for the cases of unpolarized 25|, 



26 



le target mass corrections have 



Deen studied 



27| and polarized [281, l29|, |30|, 



31 



32l | nucleon 

structure functions. As for the polarized virtual photon structure functions gJ{x,Q'^,P'^) 
and gJ{x,Q'^,P'^), TME have been studied in Ref. |33 |. 

In the present paper, we investigate the TME for the unpolarized virtual photon structure 
functions, F^{x, Q\ P^) up to the NNLO and P2(x, Q\ P^) up to the NLO, in pQCD. We use 
the framework of the OPE supplemented by the RG method. The photon matrix elements of 
the relevant traceless operators in the OPE are expressed by traceless tensors. These tensors 
contain many trace terms so that they satisfy the tracelessness conditions. The basic idea for 
computing the target mass corrections is to take account of these trace terms in the traceless 
tensors properly. There are two methods used so far for collecting all those trace terms. 
One, which was introduced by Nachtmann 25|, is to make use of Gegenbauer pa 



ynomia. 



to express the contractions between q^^ ■ ■ -q^^^ and the traceless tensors [25|, 128|, 129|, 130 1. 
This method leads to the Nachtmann moments for the operators with definite spin. The 
other, first used by Georgi and Politzer 26|], is to write traceless tensors explicitly and then 
to collect trace terms and sum them up. Through the latter approach, the moments of 
structure functions are expressed as functions of the reduced operator matrix elements and 
coefficient functions with different spins. Actually both methods give equivalent results. In 
this paper we apply the former method to study the target mass corrections to the structure 
functions P^ and P^. 

In the next section we discuss the framework for analyzing the TME based on the OPE. 
We introduce Gegenbauer polynomials to take account of the trace terms properly. In section 
3 we derive the Nachtmann moments for the structure functions using the orthogonality 
relations of Gegenbauer polynomials. In section 4, by inverting the Nachtmann moments, 
we obtain the explicit expression for F2 {x , Q"^ , P"^) (for F2{x,Q'^ , P"^)) evaluated up to the 
NNLO (up to the NLO) with TME included. In section 5 we perform the numerical analysis 
and show that target mass corrections become sizable near x^^x- The final section is devoted 
to the conclusion. 
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FIG. 2: Forward scattering of a virtual photon with momentum q and another virtual photon with 
momentum p. The Lorentz indices are denoted by fi, v, p, r. 



II. OPERATOR PRODUCT EXPANSION 

We analyze the virtual photon structure functions Q^, P^) and F2{x, Q^, P^) using 

the theoretical framework based on the OPE and the RG method. Unless otherwise stated, 



we will follow the notation of Ref. 



5j. Let us consider the forward virtual photon scattering 



amplitude for 7(g) + 7(p) 7(g) + 7(p) illustrated in FigJJl 

Tp.pr{p,q) = t [ d^xdSd''ze'^-^e'P<y''\<d\T{J,{x)JMJp{y)Jr{z)m , (2-1) 



where is the electromagnetic current. Its absorptive part is related to the structure tensor 



W^uprip, q) for the target photon with mass squared p^ = —P^ probed by the photon with 



WpupT{p,q) = -lmTp^pr{p,q) . (2.2) 

71 



Taking a spin average for the target photon, we get 
W;i,{p,q) = lY.<x^P^^^'^pr{p,q)el,){p) 



2 

A 

= -l9'^Wp,pr{p,q) = I I d'xe'^^{^{p)\Jp{x).MOMp))spln av (2.3) 

Now W^^{p,q) is expressed in terms of two independent structure functions F2{x,Q'^,P'^) 
and F2{x,Q'^,P^) without neglecting the target mass squared (see Appendix IXj) : 



(2.4) 
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where 



d^u = -g^u + I , 2.6 

p.q [p ■ q) 

and X is the Bjorken variable defined by x = Q'^/2p ■ q. 

Applying OPE for the product of two electromagnetic currents at short distance we get 

n— even 

E (l^)^M. EC'2,na^^^ -""-^ + ■■■ , (2.7) 



n=2 
n— even 



where „ and „ are the coefficient functions which contribute to the structure functions 
F2 and respectively, and 0'^^ ^" and O^"^^ ^""^ are spin-n twist-2 operators (hereafter 
we often refer to Of^' ^" as O"). The sum on i runs over the possible twist-2 operators and 
■ ■ ■ represents other terms with irrelevant coefficient functions and operators. In fact, the 
relevant O" are flavor singlet quark {ip), gluon (G), flavor nonsinglet quark {NS) and photon 
(7) operators. It is noted that the operators O" are traceless and have totally symmetric 
Lorentz indices /xi ■ ■ ■ /i„ (Acr/ii ■ ■ ■ fin~2)- 

The spin-averaged matrix elements of these operators sandwiched by the photon states 
with momentum p are expressed as 

(7(p)|Of^-^"|7(p))spinav = AUfi\P^) - trace terms} 

^ A^(/i^p2)|p,,...pM.}^^ (2.8) 

where i = ip,G, NSj'j, and y4^(/x^,P^) is the reduced photon matrix element with /x be- 
ing the renormalization point which we choose at /i^ = P^. For —p'^ = ^ A^, 
we can calculate y4^(P^) perturbatively. The {p^^ ■ ■ ■ p'^"}n denotes the totally symmet- 
ric rank-n tensor formed with the momentum p alone and satisfles the traceless condi- 
tion g^i^jip'^^ ■ ■ ■p^"}n = . Taking the spin-averaged photon matrix elements of (12.71) we 
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obtain for the photon-photon forward-scattering amphtude 



T;,(p,g) = z / dW^-(7(p)|T(J^(a;)J.(0))|7(p)) 



spin av 



n— even 



E {AV^---p^"-^}n 5^ Q,„AUP2) . (2.9) 



n^2 
n— even 



The basic idea for treating target mass corrections exactly is to take account of trace terms 
in the traceless tensors properly. We evaluate the contraction between q^^ - ■ ■ q^,, and the 
traceless tensors without neglecting any of the trace terms. The results are expressed in 
terms of Gegenbauer polynomials 25|, |28|, |29|] : 



q,,---q,AP^'---p'-}n = a-C^^\r^) 



(2.10) 



n{n — 1) 



^a"2Cfi(r^) + ^a"8C2,(r^) 



f2.111 



where 



f2.12) 



and Cn\ri)'s are Gegenbauer polynomials (see Appendix [B]). Recall that in the case of a 
nucleon target with mass M, we had = M^, rj = ip ■ q/MQ and a = —^iMQ In 
the photon case, we have p"^ = —P^ instead, and thus, replacing M with —iP, we obtain 
the expressions for a and rj in (I2.12p . The derivation of Eqs. fl2.10p and (12.110 are given in 
Appendix O 

We decompose the amplitude T'^yip^q) as 



X {p ■ qY X 



X 



(2.13) 
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then, using the results (12.10p and (12. lip , we find from Eg. (12.90 

1 



n^2 
n— even 



n(n — 1) 



2,n 5 



(2.14) 



1^7 I P q 1 

X ^ {p ■ qY X 



P 

Q 



n— even 



n— 2 
n— even 



Q / ^(^ — 1) 



X 



d22{v)-'^C':Uv)+^vC':iM Ml^, (2.15) 



.(3) 



(3) 



where we have defined 



III. NACHTMANN MOMENTS 



We derive the Nachtmann moments for the definite spin-n contributions, and Ml^. 
First we write the dispersion relations for and T2, and we denote 



F2 = -Im (p ■ q)T. 



2 ■> 



Fl = -ImT2 . 
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(3.1) 



Then, using the orthogonality relation flB4p and the integration formula fIBSP for the Gegen- 
bauer polynomials Cn\f])j we can project out and M2„. The results are as follows: 



f^UQ^P') = / dx-c-"' 

Jo x-^ 



3 + 3(n + l)r + n{n + 2)r^ 
(n + 2)(n + 3) 



a;' 



+ 



4PV (n + 3) - (n + l)e^PVQ^ 
~Q2 (n + 2)(r2 + 3) 



(3.2) 



Ml^. (3.3) 



The Nachtmann moments filn /i^n given by the weighted integrals of the structure 
functions P^ P^ and are equal to the definite spin-ra contributions, M^^ and M2„, 



respectively. The variables r and ^ are defined as 



1 



4P2x2 



2x 



2x 



Q 



2 ' 



1 + 1 + r 



(3.4) 



We see from Eq. (ll.2p that the maximal value of x is not 1 but 1/[1 + {P'^/Q'^)]. Therefore, 
the allowed ranges of r and ^ turn out to be r^i^i r < 1 and < ^ < 1, respectively, where 

Tmin = r(w) = (1 - PVQ')/{1 + PVQ') and e(Xn,ax) = 1- 
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We now outline how to derive the Nachtmann moments for the case of F2{x,Q'^,P'^) 
given in fl3.2p . Since (p- q)/x = {—PQri){—P/Q)2ri, we see that T2 in fl2.14p is expressed as 

n.-2 H 

(3.5) 



Ti — 2 

Ti— even 



By the use of orthogonahty relation of the Gegenbauer polynomials (lB4p for z/ = 3 we get 



n~2 



;i - v'f'ci^Uv) T^dv = (-^ j (n + 2)(n + 3)M2^,„ . (3.6) 

Applying the dispersion relation, we can relate the full amplitude with its absorptive part: 



1 1 

+ 



1 



Xr, 



(3.7) 



where u = 2p ■ q/Q^ = 1/x and denoting ( = (u;'Q)/(— 2P), r] = (luQ) / {—2P) we derive 



-1^ 



1 Ju^in [~^-^) J-1 [~^) ~ V 



l-r,yC^^lMdri 



(3.8) 



where we have noted C^}_2{~V) = ^n-iW) even n. From fIBSP for m = 0, z/ = 3, we get 



.(3) 



(1 - ri'f"C^:i2{v)T:idr^ 



-1 



2P / 4 



da;iy2''(a;, Q^ P2) (C^ - 1) (-1) [-C - (C' - 1)'/'] {n + 2)(n + 3) 



X 



1 



6 



-z + 



12 



(3.9) 



n + 2 (ri + 2)(n + 3) 
where we used the following relation for the hypergeometric function with C, = —p ■ q/PQ 



F{3,-2;n + 2;z) = 1 



6 



12 



,^z(z|)±Klzi)^.(3.io) 



n + 2~ ' {n + 2){n + 3y ' ~ 2{C - ly/'^ 

Setting (13. 6p equal to (13.90 and changing integration variable from u to x, we get 



dx—C^^ I 1 
x-^ 



6 



:z + 



12 



n + 2 (n + 2)(n + 3) 



F^{x,Q\P')=Ml^. 



(3.11) 



Here one should note that n is even and the following relations hold: 

Q A //•2 iNi/2 -P^ 



c 



2Pa;' 



-C - (C - 1 



2P2x 



:i-r) 



(3.12) 



and so we finally get for Q^, P^) as 



dx—^ 



n+l 



3 + 3{n + l)r + n{n + 2)r 



21 



F^ix,Q',P') = MUQ',P') , (3.13) 



(n + 2)(n + 3) 

where the left-hand side is /ig „(Q^, -P^), the Nachtmann moment, which is equal to the 
definite spin-n contribution , and this is consistent with the previous result for the case 



of nucleon target [25 



34 



35 



36| with a replacement of the variable M —iP. 
For the longitudinal structure function, we first solve for -T2 from Eqs. fl2.14l) and fl2.15p 



and we get (see Appendix iBl) 



x 



■ rp"f 



E 

n— even 



P 
Q 



n=2 
n— even 



Q J n{n — 1) 



n 



n 



Mir. 



Then using the recursion relation flBlOp for the case v = 1: 



(3.14) 



Cl^lM , (3.15) 



(n + l)(n + 2) n{n + 2) " n{n + 1) 

(3) 

and the orthogonality relation of the Gegenbauer polynomials C„ s we can derive the re- 
cursive relations for the sequence M2„'s which can be solved as (See Appendix [D|) : 



Jo 



dx—i 

X'^ 



71+1 



'Flix. Q\ P^) + l^ (n + 3)-(n+l)eJ/g^ ^,^^^ 



(n + 2)(n + 3) 



(3.16) 



where the left-hand side integral is the Nachtmann moment /i^^ for the longitudinal part. 



This coincides with the result obtained in 



34( 1 after the replacement mentioned above. 



IV. INVERTING THE NACHTMANN MOMENTS 

We now invert the Nachtmann moments to express the structure functions and 
explicitly as functions of x, and P^ . We first consider F^ - By changing the integration 
variable from x to ^, we can rewrite the Nachtmann moments, given in Eqs. (13.21) and (13.31) . 
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as follows: 



ML 



^0 



X 



3 3(n+l) n /l-n^'^^-^ 

(n + 2)(n + 3) ^ (n + 2)(n + 3) 1 + + ^ + 31,1 + ^^2 



F2'^(a;,g^O (4.1) 



^7. 



JO 



X 



F2(x,Q^P^) + 4«: 



n + 1 



(1 + K^'^)'^ \n + 2 ^ (n + 2)(n + 3) 
where we have made use of the following relations: 



(4.2) 



X 



l-ne 



dx 



with K — P'^/Q'^. Now we define 

1-^:^2 



1 + <2 ' (1 + nef ' 



(4.3) 



f:j{x,q\p'') , B{0^{i-Ke){i + Ke)F2{x,Q\p'^ 



1 + 

Then the above two moments are written as 



(4.4) 



ML 



n-2 



A(0 2KeA{0 , «:r^(0 



n(n — 1) Jo \n{n — 1) n{n + 2) (n + 2)(n + 3) 

'0 



M7„ 



(4.5) 
(4.6) 



n + 1 Jo \n + l ' (n + l)(n + 2) (n + 2)(n + 3) 

The boundary conditions for A{^) and B{^) are = 1) = = 1) = 0, since 
Fiixmax^Q'^^P'^) = -F2(^max, <5^ -P^) = and C(xinax) = 1- Now introducing the following 
four functions, 



s_3(0 = f dc 



(4.7) 
(4.8) 



and by partial integration we find that the above two moments are written as 



n{n — 1) 



Ml, 



= ^""'(1 - + - ^^-2(0 - <^o(e)}, (4.9) 

I = ^'c?Cr[5-3(0+4«:{(l + «:aA_i(0-e^-2(0-<^(0}]- (4.10) 
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Inverting the moments we get 



-| pc+ioo f 



ML 



n{n — 1) 

:i - + Ke)A.,{o - e^-2(o - <Aoio}, (4.11) 



1 /"C+ioo I' /\/f7 

1 / . ^^^L,n 



S{0 = TT- dni 



4k 



where in Eq. (14. 12^ we have used the resuh of Eq. (I4.1ip . 
We further introduce the following functions, 



H{0 ^ 


dG{0 
d^ 


2ni J 


'c-ioo ^ 


(4.13) 


m = 


dH{i) 
di 




fc+ioo 

/ dnC^'-^Ml^ 

1 c—ioo 


(4.14) 




^ di 


1 

" 2TTi 


fc+ioo 

/ dnC'-'Ml^ . 

J c—ioo 


(4.15) 



Differentiating both sides of Eqs. f l4.1Hl4?T2l) with respect to ^, we get the relations between 

A_i(0, (A_2(0 + 't^o(O), ^(0, B{0 and G(0, H{0, F^iO and ^^(0. Now solving for 
A[^) and B{^) and, then recalling Eq. (14.41) . we obtain 



2 4 

rf-" rr"-^ 



F2^(x, P') = ^FiO - Qf^^HiO + l2K'^GiO , (4.16) 



^3 



F2(x, Q^ P2) = ^F^(0 - 4«:^i/(0 + SK^^GiO . (4.17) 



Eqs. (l4.16p . (I4.17P are the final formulas for the photon structure functions and F2 
when target mass effects are taken into account. They can also be derived from the method of 



Georgi and Politzer 26|]. Once and in Eq.(2.16), are given, then we can calculate 
the four "profile" functions G(0, ^^(0. ^(0 and Fl{0 through Eqs. fl4. 11114.151) . and by 
using Eqs.f l4.16ll4.17p we can predict whole structure functions with target mass corrections. 

Note that in the above expressions, the Q^- as well as P^- dependence of F{^), G{^), 
iJ(^) and -Fl(0 are implicit, since they are given by and M2„ which depend on as 
well as on P^. If we take the k —>■ limit, the above expression reduces to that for the F2 
and F2 without TME. And in the absence of TME we have 



1 



ML 



/'dxx"-2F2^(x,Q^p2), M2^^= f dxx''-^F2{x,Q\P^) . (4.18) 
Jo ' Jo 
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From our previous QCD calculation 21j of M]^ and M2„ we already know the three 
functions, F{C,), G{^) and H{C,) to NNLO and to NLO, so we can evaluate the photon 

structure functions with TME to the same accuracies as in the case we neglect TME. 

V. NUMERICAL ANALYSIS 

In this section we perform a numerical analysis for the structure functions and 
when TME are included. We first compute the four profile functions F, G, H and Fl which 
are given in Eg. (HIT]) and Eqs.dHSD-dHSD- We use the QCD results for M^^ and M2_, 
which have been calculated up to the NNLO and the NLO in QCD, respectively, in Ref. 21|. 



Indeed, the expressions of and M2„ are given in the right-hand sides of Eq.(2.29) and 
Eq.(6.3) of Ref. {21!]. In Fig. 3 we have plotted the functions F, G, H and Fl as functions 
of X for the case of = SOGeV^ and = iGeV^ with x^ax = 0.968. We take A = 0.2 
GeV for the QCD parameter and Uf = 4 for the number of active quark flavors throughout 
our numerical analysis. Note that we have multiplied each function by a suitable power of 
X to accommodate four functions in a single graph. The Bjorken variable x ranges from 

to ^^max- 

Now inserting the functions F, G, H and Fl into Eqs. (14.161) and (14.171) . we obtain the 
graphs of -^^(a;, Q^, P^) and F2{x, Q^, P^) as functions of x, which are shown in Fig. 4 and 
Fig. 5, respectively. We observe that TME become sizable at larger x region. While TME 
enhances F2 at larger x, it reduces F2- In fact, F2 becomes maximum at x very close to the 
maximal value of x, Xmax (1) for the case with (without) TME. The target mass correction is 
of order 10 % when compared at the maximal values for F2 . In the case of F^, the maximal 
value is attained in the middle x, where the TME reduces the F2 about 5 %. 



One should note that P(0 is dominant at larger x region, and determines the leading 
behavior of F2 ■ And the factor x^ jr^ in front of F2 shows the deviation upwards from x^ 
as X approaches Xmax- On the other hand, the expression for F^-, Eq. (l4.17l) . possesses no 
dependence upon the function P(0, in contrast to the case of F^- This is a reason why P^ 
becomes maximum in the middle x region, as seen from Fig.3. 
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FIG. 3: The four functions F, G, H and Fl as functions of x. = soGeV^ and = iQeV^. 
a^max = 0.968. 



Q. 

CM 
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FIG. 4: F^{x, Q^,P'^) as a function of x for = SOGeV^ and = IGeV^ with Xmax = 0.968. 



Now let us compare our theoretical prediction for the virtual photon structure functions 
with the existing experimental data. In Fig. 6 and Fig. 7, we have plotted the experimental 



data from PLUTO Collaboration 



371] and also those from L3 Collaboration 38|] on the so- 



called "effective photon structure function" defined as F^^ = F2 + §-^2; together with the 
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n. = 4, Q ^SOeV^, P'^=0.35GeV'^, x 



NNLO + no TME 
NNLO + TME 
BOX(LO + NLO) 
PLUTO 



Q. 




FIG. 6: NNLO predictions with and without TME for the effective photon structure function: 



F^r. = + |f7 for Q2 = 5GeV2 and = O.SSGeV^ with 



0.93. The experimental data 



are from the PLUTO group [37]. The Box diagram prediction to the NLO order is also shown 
(blue short-dotted line). 



15 



1.5 



eg 

o 



0.5 



n, = 4, Q^=120GeV^,| P^=3.7GeV^, x^ax=0.97 



NNLO + no TME 
NNLO + TME 
BOX(LO + NLO) 
L3 




0.1 0.2 0.3 0.4 



0.5 
X 



0.6 0.7 



FIG. 7: NNLO predictions with and without TME and Box prediction (NLO) on Fj^ = +1^2 
for Q2 = i20GeV2 and = s.TGeV^ with x^ax = 0.970. The data are from the L3 Group [38]. 



theoretical predictions. The effective structure function is proportional to arr + o'LT + o'TL + 
aLL, where aab (a, b = T,L; T = transverse and L = longitudinal) is the total cross section 
with the helicity state (a) of the probe photon and helicity (b) of the target photon. This 



combination of and F2 is obtained in the limit of P'^/Q'^ <^ 1 [lO|, l39|, |40 |. 

In the above experiments we have = 5 (120) GeV^ and = 0.35 (3.7) GeV^ with 
a^max = 0.93 (0.97), for PLUTO (L3) data. Here we note that for the PLUTO data, P^ < 
is satisfied since P'^/Q'^ ~ 0.07, but is not much larger than i.e. A^/P^ ^ 0.114. 
For L3 data, both hierarchical conditions are satisfied. Although the experimental error 
bars are rather large, the data are considered to be roughly consistent with the theoretical 
expectations, except for the larger x region in the case of L3 data. Note that the TME for 
the in this kinematical region is almost negligible for both cases. This could be explained 
as a consequence of the cancellation of the TME between F2 and P2, as discussed above. 
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VI. CONCLUSIONS 



We have investigated the target mass corrections for the unpolarized virtual photon 
structure functions F^{x, P'^) and Fl{x, P^) to the NNLO in perturbative QCD. 

In contrast to the case of the nucleon target, the virtual photon target provides us with 
the unique testing ground for the perturbatively calculable target mass effects. Taking into 
account the trace terms in the operator matrix elements by using the expansion in terms of 
orthogonal Gegenbauer polynomials we get the Nachtmann moments. These moments were 
then inverted to derive exphcit expressions for F2{x, Q^, P^) and Fl{x, Q^, P^) in terms of 
the four profile functions which are calculable to NNLO for the first three functions P, G, 
H and to NLO for the last one P^. The TME becomes sizable at larger x region, and it 
enlarges F2 near Xmax and reduces P^ in the region of x larger than the middle point. When 
we go to higher values of P^, e.g. P^ = 3GeV^ for = 30GeV^, it has turned out that 
the F2 blows up as x approaches Xmax- So some prescription like resummation of large logs 
would be needed to avoid such difficulties. 

We have carried out the confrontation of our theoretical predictions with the existing 
experimental data on the effective photon structure function F^q from PLUTO and also 
those from L3 Collaboration. Roughly speaking we find the rather good agreement between 
theory and experiments. However, it turned out that TME looks almost negligible for P^g, 
which is the combination of F2 and P^ and exhibits a cancellation of TME between them. In 
the present analysis, we have treated the active flavors as massless quarks, and ignored the 
mass effects of the heavy flavors, which should remain as a future subject. We should also 
investigate the power corrections (P^/Q^)*^ (/c = 1, 2, • • • ) due to the higher- twist effects. 

We expect the future experiments would provides us with more accurate data for the 
double-tag two-photon processes in e+e~ collisions. 
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APPENDIX A: PHOTON STRUCTURE FUNCTIONS 



Averaging the structure tensor W^j^upriPi l) given in Eg. (12. 2 p over the target polarization, 



we get 



41 



R 



where 



R 



k 



X 



p-q 



(Al) 

(A2) 
(A3) 



with X = [p ■ qY — q^p^. In the above equation the first index (a = T, L) of the invariant 
functions Wab refers to the probe photon and the second one {h = T, L) to the target photon, 
and the subscripts T and L denote the transverse and longitudinal photon, respectively. 
We define the unpolarized photon structure functions and F]^ as. 



^2 



{p_qY 
X 



Wtt - - Wtl 



+ 



Wlt - -Wll 



Wlt - i^Wll 



Another structure function is often used, which is defined as j43| 



F7 = Wtt - -Wtl- 



Then we get a well-known relation 



Fl = -xF^ + 



X 



-xF7 + 1 - 



^2 



(A4) 
(A5) 

(A6) 

(AT) 



(p-g)2 2 

Since R^u and ki^kiy are expressed in terms of e^y and d^y, which are given in Eqs. (12.51) and 
(EID, as 

{p-qf 



R 



p2g2 

^t^. ■ X 



kifxkiy 



X 



2 



-d 



we find that W^y{p, q) in ( lAll) is rewritten as 



WUp.q) 



1 „^ p^q^ 1 



-Fl'rj-^^-F:i\^d,y-F, 
X [p ■ qY X 



X 



(A8) 
(A9) 

(AlO) 
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APPENDIX B: GEGENBAUER POLYNOMIALS 



The Gegenbauer polynomials are defined through the generating function given by 441.145 1 



(l-2vt + tY' = J2ci^\vr. (Bl) 



n=0 



In terms of hypergeometric functions F{a,(3,'y; z), Cn \ri) is expressed as 
. ), , 2"r(n + z/) / n 1-n ^ 1 



n/2 



r(f)^ j!(n-2i)! 

For example, we have 



^ j! (n-2j)! 



1. Orthogonality relations 



The orthogonality relation reads 



2^^ (n + z/)n![r(z/)]2' 



In addition we have the following formula for the integral to project out the contributions 
of definite spin from the dispersion relations, 

dr, r^-ii - vr-^c^:\v)^ = ^.ne - 1)^ [c - (c^ - ly^T" 

T(n + 2iy) ( _^ + (^2 _ i)i/2\ 

In fact the factor — [C^ — l)V2j"+'^ gives ■, where ^ is the so-called ^-scaling 

variable given in Eg. (13.41) . 



2. Recursion relations 



The recursion relations for Gegenbauer polynomials read 

nC^:^ iv) = 2u[r^Ct\'^ (r^) - C^V^ iv)] , (B6) 

(n + 2u)C^:^ iv) = MC^r'^ iv) - vCt\'^ iv)] , (B7) 

(n + 2)C(%(r^) = 2{n + u + l)vd:lM - {n + 2v)&:\r,) . (B8) 
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We get from and fiBTll 



V 



n + u 

V 



v + l 



n + V n + u + 



(n + z/)2 - 1 



n + u n + u — 1 



(B9) 



(BIO) 



Now we derive the second line of Eq. fl3.14p . Choosing u = 2 and n ^ (n — 2) in (lB9p . we 



find 



n 



(3) 



.(3) 



Next choosing z/ = 3 and n ^ (n — 4) in (]B8I) . we get 

^vC^alM = -\{n- 2)d:iM + {n + 2)C'^1m 



(Bll) 



(B12) 



Thus we obtain 

.(2) 



c):Uri) - + ^vc!^1m] - ^c^^Uv) = —ci^Uv) + -ci^lM (B13) 

71/ 71/ 



.(3) 



.(3) 



APPENDIX C: CONTRACTION FORMULAS 



Here we derive Eqs.f l2.10l) and (12.111) . The most general rank-ra symmetric and traceless 
tensor, '''P^"}n, that can be formed with the momentum p alone, is expressed as 



follows 



{pM^ . . = ^ ^—L^JL_lIlg . . . gf^ (p2), 



j=0 



23 n\ 



(CI) 



where g ■ ■ ■ g stands for a product of j metric tensors (yi^'^fc with 2j indices chosen among 

j 

/ii, ■ ■ ■ , /in in all possible ways. Then we easily see that the contraction of {p^^ ■ ■ ■ p^"|„ with 
Q/j-i ' " QiJ-n is expressed in terms of Gegenbauer polynomial Cn\rj) given in (1B3P as 25, 3], 



QtJ.1 ' ' ' QiJ., 



^{pMi...pM"} = a"C«(r^) 



(C2) 



which is Eq. (l2.10l) . Here we have put 



a = -l^PQ^ V = 



PQ 



(C3) 
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Next we differentiate both sides of Eq. (lC2p twice with respect to qa and g/3. The left-hand 
side becomes 



d d 



{q,. ■■■q,^ {P''' ■■■P"" }n) =n{n-l)q^,--- q,^_, {p^Y' " " " P^^' }n ■ (C4) 



dqa dqp 

Differentiation of aP'Cn^ {rj) with respect to q^ gives 



(C5) 



where we have used the following formulas 



da f—q" 
a 



dn rf 



2^^CtV^(^) , (C6) 



and, at the last line, the recursion relation flB6l) . Further, we differentiate the both sides of 
Eq. llCSP with respect to qp . Again using the formulas in (]C6|) and the recursion relation 
dHlD, we get 

+pVa"-= ■ 2Cf + p!9!±9>!„«-i . icflM (C7) 
Thus, from Eqs. flC4p and (\C7\} . we obtain 



n{n — 1) 



(3) , p"g'^ + gV„„-i 



(C8) 



which is Eq.(^JB>- 



APPENDIX D: NACHTMANN MOMENTS OF 

By applying the orthogonality relation ( 1B4I) for z/ = 3 to the longitudinal amplitude ( ]3.14p 
with the help of (13.151) we derive the following recursive relation for „'s: 
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where 



+Ak 



3 + 3{n + 3)r + {n + 2){n + A)r' 
(n + 4)(n + 5) 

2 (n + l)(n + 2) 



(n + 3)(n + 4) 



2 ^*^^2,n+4 



The above recursive equation (IDll) can be solved as an infinite series: 



[n 



- , (/ + l)(n + 2 + /) ^ 



i=0 



(n + 2/ + l)(n + 2/ + 2)' 



Introducing a variable t defined by ^ (P/ Q) =t we have 



1 -t^ 



dx 



1-f 



and then in terms of t we can sum up the above infinite series and find 

1-e 



_x^ " (n + 2)(n + 3) 

the right-hand side of which turns out to be the Nachtmann moments (13. 3p . /i^n- 



(D2) 



(D3) 



(D4) 



(D5) 



[1] T.F. Walsh, Phys. Lett. 36B, 121 (1971); 

S.J. Brodsky, T. Kinoshita and H. Terazawa, Phys. Rev. Lett. 27, 280 (1971). 
[2] T.F. Walsh and P.M. Zerwas, Phys. Lett. 44B, 195 (1973); 

R.L. Kingsley, Nucl. Phys. 60, 45 (1973). 
[3] N. Christ, B. Hasslacher and A.H. Mueller, Phys. Rev. D6, 3543 (1972). 
[4] E. Witten, Nucl. Phys. B120, 189 (1977). 
[5] W.A. Bardeen and A.J. Buras, Phys. Rev. D20, 166 (1979); 

Phys. Rev. D21, 2041 (E) (1980). 
[6] R.J. DeWitt, L.M. Jones, J.D. Sullivan, D.E. Willen and H.W. Wyld, Jr., Phys. Rev. D19, 

2046 (1979); Phys. Rev. D20, 1751(E) (1979). 
[7] M. Gliick and E. Reya, Phys. Rev. D28, 2749 (1983). 

[8] K. Sasaki, Phys. Rev. D22, 2143 (1980); Prog. Theor. Phys. Suppl. 77, 197 (1983). 
[9] M. Stratmann and W. Vogelsang, Phys. Lett. B386, 370 (1996). 
[10] M. Gliick, E. Reya and C. Sieg, Phys. Lett. B503, 285 (2001); 
Eur. Phys. J. C20, 271 (2001). 



22 



[11] T. Uematsu and T. F. Walsh, Phys. Lett. lOlB, 263 (1981). 
[12] T. Uematsu and T. F. Walsh, Nucl. Phys. B199, 93 (1982). 
[13] G. Rossi, Phys. Rev. D29, 852 (1984). 

[14] M. Drees and R. M. Godbole, Phys. Rev. D50, 3124 (1994). 

[15] M. Gliick, E. Reya and M. Stratmann, Phys. Rev. D51, 3220 (1995); Phys. Rev. D54, 5515 
(1996). 

[16] M. Fontannaz, Eur. Phys. C 20, 297 (2004). 

[17] K. Sasaki and T. Uematsu, Phys. Rev. D59, 114011 (1999). 

[18] K. Sasaki and T. Uematsu, Phys. Lett. B473, 309 (2000); Eur. Phys. J. C20, 283 (2001). 

[19] K. Sasaki, T. Ueda and T. Uematsu, Phys. Rev. D73, 094024 (2006). 

[20] M. Krawczyk, AIP Conf. Proc. No. 571 (AIP, New York, 2001) and references therein; 

M. Krawczyk, A. Zembrzuski and M. Staszel, Phys. Rep. 345, 265 (2001); 

R. Nisius, Phys. Rep. 332, 165 (2001); 'hep^ex/0110078 



M. Klasen, Rev. Mod. Phys. 74, 1221 (2002); 

I. Schienbein, Ann. Phys. 301, 128 (2002); 

R. M. Godbole, Nucl. Phys. B (Proc. Suppl.) 126, 414 (2004). 
[21] T. Ueda, K. Sasaki and T. Uematsu, Phys. Rev. D75, 114009 (2007). 
[22] S. Moch, J.A.M. Vermaseren and A. Vogt, Nucl. Phys. B688, 101 (2004). 
[23] A. Vogt, S. Moch and J.A.M. Vermaseren, Nucl. Phys. B691, 129 (2004). 



[24] A. Vogt, S. Moch and J.A.M. Vermaseren, Acta Phys. Polon B37, 683 (2006); [hep^ph/0511112 

[25] O. Nachtmann, Nucl. Phys. B63, 237 (1973); B78, 455 (1974). 

[26] H. Georgi and H. Politzer, Phys. Rev. D14, 1829 (1976). 

[27] A. De Rujula, H. Georgi and H. Politzer, Ann. of Phys. 103, 315 (1977); Phys. Rev. D15, 
2495 (1977). 

[28] S. Wandzura, Nucl. Phys. B122, 412 (1977). 

[29] S. Matsuda and T. Uematsu, Nucl. Phys. B168, 181 (1980). 

[30] H. Kawamura and T. Uematsu, Phys. Lett. B343, 346 (1995). 

[31] A. Piccione and G. Ridolfi, Nucl. Phys. B513, 301 (1998). 

[32] J. Bliimlein and A. Tkabladze, Nucl. Phys. B553, 427 (1999). 

[33] H. Baba, K. Sasaki and T. Uematsu, Phys. Rev. D65, 114018 (2002). 

[34] S. Simula, Phys. Lett. B574, 189 (2003). 

23 



[35] F. M. Steffens and W. Melnitchouk, Phys. Rev. C73, 055202 (2006). 

[36] For a recent review, see for example, I. Schienbein et al., larXiv: 0709 .17751 [hep-ph] . 

[37] PLUTO Collaboration, Ch. Berger et al., Phys. Lett. B142, 119 (1984). 

[38] L3 Collaboration, M. Acciarri et al., Phys. Lett. B483, 373 (2000). 

[39] R. Nisius, Phys. Rep. 332, 165 (2001). 

[40] M. Gliick, E. Reya and I. Schienbein, Phys. Rev. D63, 074008 (2001). 

[41] V. M. Budnev, V. L. Chernyak and I. F. Ginzburg, Nucl. Phys. B34, 470 (1971). 

[42] K. Sasaki, J. Soffer and T. Uematsu, Phys. Rev. D66, 034014 (2002). 

[43] Ch. Berger and W. Wagner, Phys. Rep. 146, 1 (1987). 

[44] A. Erdelyi et al.. Higher Transcendental Functions, Bateman Manuscript Project, Vol.1 

(McGraw-Hih, New York, 1953). 
[45] I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series, and Products (Academic Press, 

New York, 1965). 



24 



